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Abstract 

The study of the convergence of power series expansions of energy 
eigenvalues for anharmonic oscillators in quantum mechanics differs from 
general understanding, in the case of quasi-exactly solvable potentials. 
They provide examples of expansions with finite radius and suggest tech¬ 
niques useful to analyze more generic potentials. 
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1 Introduction 


Let me summarize some of the most relevant features of perturbative expansions 
for eigenvalues in one-dimensional quantum mechanics. Let us consider the 
Schrodinger equation 


-— + - + \V{x) 


ipk{x) = Ekipkix) 


( 1 . 1 ) 


where the harmonic oscillator is perturbed by a higher order even polynomial 
XV(x) = Xp{x^) and k is the number of zeros of the wave function . 

In well known papers C.Bender and T.T. Wu 0 and B.Simon [|| studied the 
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quartic anharmonic oscillator p{x^) = . The former authors Q| evaluated a 

large number of terms of the perturbative series 

OO 

T;fc(A) = (fc + l/2) + ^A"4”^ (1.2) 

n—1 

of the lowest energy eigenvalues Ek{\), k = 0,1,.. and discussed the occurrence 
of singularities in the complex A plane. It was found that for any energy level 
EkiX) there exist infinitely many points A , which are extrema of square root 
branch singularities , the extrema accumulate at the origin in the complex A 
plane and each value A correspond to the crossing between pairs of energy levels 

E„iX)=EM (1.3) 


These features prevent a non vanishing radius of convergence for the perturba¬ 
tive expansion of any energy eigenvalue ( 0 )- The second author confirmed 
these results, by using Hilbert space methods. Analogous occurrence of infinitely 
many singularities with accumulation point at the origin, for any energy level, 
was found if the quartic monomial was replaced by a higher order mono¬ 
mial Ax^", n > 2 [^. T.Banks and C.Bender i also studied the anharmonic 
oscillator with a general polynomial potential (of even parity) 


9 ^ 



(1.4) 


The k — th energy level has a perturbative expansion 

- OO 

Eu{g) = k + - + Y.^n\aA-)9'^ 

n—\ 


They found that the large order behaviour of the coefficients c^\a, b ,..) is given 

by 


Cn\a,b ,..) 

ci'=^(0,0,..) 


ga/(7V-l) 


l + 0(-) 

n 


(1.5) 


Then the singularities of Ek{g) closest to the origin in the complex g plane 
are controlled by the highest order monomial g{x‘^/2)^ whereas the next or¬ 
der monomial ga{x‘^/2)^~^ only results in a constant factor and the next order 
monomial gb{x'^/2)^~‘^ affects the corrections of order 0(l/n) with respect to 
the previous result. 

For some decades it was believed that any formal Taylor expansion of energy 
eigenvalues of an anharmonic potential, with any polynomial perturbation (of 
degree higher than quadratic) would have a vanishing radius of convergence. It 
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was recently found that for the class of models known as quasi exactly solv¬ 
able potentials, a number of energy levels have a perturbative expansion with 
finite radius of convergence Q . The singularities of these energy levels still cor¬ 
respond to level crossing, yet these are a finite number. Quantum mechanics 
being a (0-1-1) dimensional quantum field theory, it would be exciting to find 
similar convergence in higher dimensional models of quantum field theory. 
Further references to extensive investigations on the divergence of the pertur¬ 
bative expansion in quantum mechanics and in quantum field theory may be 
found in ^ and j^. 

It is clear that quasi-exactly solvable models have convergent perturbative ex¬ 
pansions for a finite number of energy eigenvalues because those eigenvalues 
are decoupled from the rest of the spectrum. Yet this property is so peculiar, 
that it is interesting to have a pattern of the radius of convergence as function 
of the parameters. This is evaluated in sect.2, in an algebraic exact fashion, 
for a simple sequence of potentials. It also seems that quasi-exactly solvable 
potentials provide efficient tools to investigate the possibility of convergence for 
generic potentials, when these are summed in a fashion similar to quasi-exactly 
solvable models. This analysis is presented but not completed in sect.3. 


2 Quasi-exactly solvable potentials. 


The generic conclusion of divergence of perturbative expansion does not hold 
in the case of quasi-exactly solvable potentials. This is clearly stated in the 
book m . In this section it will be exhibited by the evaluation of the radius of 
convergence in a sequence of cases. 

The simplest class of quasi-exactly solvable potentials corresponds to the one¬ 
dimensional quantum sextic oscillator model with hamiltonian 


^2 

Hm = ~'q~2 ~ a(4M -I- 3)]x'^ + 2abx'^ -I- a?x^ 


( 2 . 1 ) 


where a is positive, b is real, M is a non-negative integer (M = 0,1,..) For sake 
of a simpler exposition, let us choose 6=1 (which is a generic value). It can be 
shown that the eigenvalue equation 


- 7 ^ + [1 — a(4M -I- 3)]a;^ -I- 2ax‘^ + a^x^] 
ox^ 


ipkix) = Eki’kix) (2.2) 


where ipk(x) is square integrable, has the lowest part of the spectrum corre¬ 
sponding to the even wave functions, which may be computed in closed form in 
algebraic way. That is, the first M+1 even wave functions ip 2 k{x'^), k = 0,1, ..M 
are 


i’2k{x^) 



(2.3) 
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(2.4) 


n—Q 


{2n)\ 


where ^i, •.. ,^m are real numbers satisfying the system of M algebraic equa¬ 
tions 


M 

E 


^ + E - 1 - 2aC* = 0 


4^i 


(2.5) 


Each of the M + 1 solutions of the set {^i} is characterized by having k of the 
numbers positive and the remaining M — k being negative. It provides one 
of the M + 1 computable energy eigenvalues: 


M 


E 2 k — {4:M -|- 1) -|- 8a 


( 2 . 6 ) 


and the ground state corresponds to the solution where all the numbers {^i} 
are negative. The eq.(2.5) has the familiar form of a saddle point equation 
for random matrix models and the promising relations between quasi exactly 
solvable models and random matrix models are just beginning to be explored 

0 §• 


The system (2.5) and eq.(2.6) imply that for hxed non negative integer M, the 
M + 1 eigenvalues of the lowest even wave functions are the roots of a polynomial 
equation of order M + 1 which may be obtained by techniques of symmetric 
functions. It is however easier to obtain them by inserting the ansatz ( [2.3| ) in the 
eigenvalue equation (2.2) bypassing the evaluation of the set {^i}. For example, 
the five polynomial equations which correspond to the values M = 0,1, ..,4 are 

E-1 = 0 

- 6E -P 5 - 8a = 0 

- 1S>E^ + (59 - 64a)E -f (192a - 45) = 0 

E^ - 28E^ + (254 - 2400)^^ - (812 - 2592a)E + (585 - 4656a -f 2880a2) = 0 


E^ - 45^"* -h 10(73 - 64a)E^ - 114(45 - 128a)E^ + 

QQ4iT 

128(^— - 709a -f 368a^)E - 128(-— - 984a -f 17760^) = 0 


128 


128 


(2.7) 


The second ansatz for the wave function (2.4) is very useful to derive the wave 
function because one immediately Hnds a three term recursion relation for the 
coefficients Pn{E) 


{E-4k- l)Pk = Pk+i - 8ak{2k - l)(fc - 1 - M)Pk-i 


( 2 . 8 ) 


with 


Pn = 1 


Pi = E-1 


(2.9) 
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The coefficient Pn{E) is then a polynomial in E of order n. The condition that 
Pm+i{E) = 0 leads to the algebraic equation for E of degree M + 1 (the lowest 
ones being eq. (EH) ). This condition and the recursion relation eq.( ^.^ ) imply 
that all Pk{E) with k > M vanish. The finite set of non-vanishing polynomials 
PniE) is a set of weakly-orthogonal polynomials, recently discussed by several 
authors §,0,0- The papers 0| , 0| show that a set of weakly-orthogonal 
polynomials occur in any quasi-exactly solvable model and they exhibit the 
discrete weight function uj{E) = ~ ^k) for which 


Pn{E)Pm{E)uj{E)dE = hnS„ 


( 2 . 10 ) 


The roots of the polynomial eqs. (|2.7| ), Pm+i = 0, M = 0,1,..., determine the 
even lowest energy eigenvalues Eo{a), £^ 2 ( 0 ), ..,£ 2 M(a)- The singularities of 
the function E 2 k (ci) closest to the origin, in the complex plane of the variable a 
provide the radius of convergence of the perturbative expansions 


£2fc(a) = l + 4fc + ^dWa" (2.11) 

n=l 


The singularities only occur for the values a such that £ 2 / 0 ( 0 ) is a multiple root 
and may be found by examining the solution of the system 

£M+i(a,£) = 0 = ^£M+i(a,£) (2.12) 

For instance, for M = 2, 

P^{E) = £3 - 15£2 -t (59 - 64a)£ -k (192a - 45) = 0 (2.13) 


defines, in closed form, the three energy levels £o(a), £ 2 ( 0 ), £ 4 ( 0 ), and their 
perturbative expansions may be easily evaluated at arbitrary order 


£o(a) = 1 — 4a — 2a^ -I- 4a^ -I- -a^ 


- lla'^ 


39 


57 7 2235 

—a - 

2 32 


a — 


563 


21809 

64 


0 ^*" + ... (2.14) 


E 2 ia) =5-8a + 32a^ - 160a^ -k 1024a^ - 75520*^ -k 59904a® - 497664a^ -t 

-h4276224a® - 37697536a® -t 3390423040^° -k .(2.15) 


2049 


£ 4 (a) = 9 -k 12a - 300^ -k 156a°- —a^ + 7563a® - 


239655 


995271 7 136836933 „ 

+ -+ 


603161139 q 21698729265 

a -rr-a -|- ... 


16 


64 


(2.16) 
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The singularities of Eo{a) , eq.( p.l3 ), occur for the three values a 

69 


= [ 3(11 + 64 ^ 3 ) 1 / 3 - 7 - 


(11 + 6473 ) 1/3 


]/64 ~ -0.09445127 


-0.116837+ 0.389587i 


(2.17) 


where the above pair of complex conjugate values are the roots of 

_2 7406 + 1587p + (33 - 19273)p2_ 

a + -- a 


33856 

p = (11 + 6473 ) 1/3 


69 + 7p - 37 


= 0 


(2.18) 


It is easy to check that the closest, real negative value of a , ( ^.17] ), corresponds 
to the radius of convergence of the perturbative expansions of both the levels 
E 2 {a) and Ei{a), (2.15) and ( ^.16 ), and may be interpreted as the value of 
a corresponding to the crossing E 2 {a) = Eiia) , while the couple of complex 
conjugate values ( ^.Ig ) correspond to the radius of convergence of the per¬ 
turbative expansion of the ground level ifo(a) , ( p.l4 ). It is remarkable how 
easily this situation generalizes for all integer values of M. The M + 1 roots 
of the polynomial equation Pm{E) = 0 d efine the M + 1 energy eigenvalues 
£'o(a), E 2 {a), ...E 2 m{o)', the system ( ^.12 ) leads to a polynomial equation in 
the variable a of degree M{M + l)/2. Its M{M + I)/2 roots in the complex a 
plane are in one-to-one correspondence with the possible level crossing among 
pairs of eigenvalues 


E2ria) = E2sia) 


(2.19) 


All these singular values were examined, for the cases M = 1 up to M = 7, 
beginning with the values a with the smallest modulus. For any M examined, 
the first value a occurs on the real negative axis in the a complex plane, it cor¬ 
responds to crossing of the two highest eigenvalues considered (for M = 7 it is 
Ei 2 (a) = Ei 4 (a)) thus providing the radius of convergence of their perturbative 
expansions ( 2.11 ). Singular values a with larger modulus describe level crossing 
between pairs of intermediate levels. Only the M values a with larger modulus 
describe level crossings of the ground state level with the other M levels. The 
radius of convergence of the perturbative expansion of the ground level E^ (a) is 
determined by a couple of complex conjugate values of a with the smallest mod¬ 
ulus in this last group of M values a. This analysis is confirmed by the study 
of the period of oscillations of the coefficients in the perturbative expansion of 
Eo{a). 


Fig.l shows a sequence of points a in the complex a plane. The point a = 
— 1/2 corresponds to the only singularity of the ground state eigenvalue ifo(a) = 
3 — 271 + 2a of the potential (2.2) with M = 1. Moving clockwise in the upper 
half plane, there are points corresponding to the singularity of Eo{a, M) closest 
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Leading Singularity of the Ground Energy 
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Figure 1: 

to the origin, for the cases M = 2,3, .., 6 . One sees that the radius tm of 
convergence decreases as M increases. 

3 Softly broken quasi-exact solvability 

The analysis of quasi-exactly solvable models is useful for more general polyno¬ 
mial potentials, where the couplings do not obey the constraints of quasi-exact 
solvability. Let us consider the hamiltonian with a generic sextic potential 


H =+ V{x) =+ ax^ + fix'^ (3.1) 
ox^ 

with 7 > 0, a and /? real. By choosing 

= (3.2) 


the hamiltonian (0 reproduces the quasi-exactly solvable model (1^), but 
now M is real, rather than a non negative integer. In this Section, I study the 
perturbation theory for the eigenvalues of this hamiltonian, by keeping b and 
M fixed, in series of powers of the coupling a. This i s ju st a specific way to do 
perturbation theory for the generic sextic potential (3.1). For simplicity let us 
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fix again b = 1, which is a value with no special meaning, and let us choose the 
(formal) ansatz for the even parity wave function 

(j){a, x^) e~~ ~ 

V (3.3) 

n (2^)- 

n—0 ' ^ 

The eigenvalue eq. €3) implies the recursion relations 

Pn{E,a,M) = (i? - 4n + 3)Pn-i{E,a, M) - A„P„_ 2 (i?, a, M) 

A„ = 8a(n-l)(2n-3)(M + 2-n) (3.4) 

Since M is not a non-negative integer, the number of non-vanishing polynomials 
Pn{E, a, M) is infinite. However A„+i > 0 only for n < M + 1 so the system of 
infinite polynomials Pn{E, a, M) is not an orthogonal systems with respect to 
a non-negative Stieltjes measure [^. Still the recursion relations (|3.4D, are a 
powerful tool for the perturbative analysis of the generic sextic potential &■ 
The first few polynomials Pn{E,a, M) are 


= 


Po{E,a,M) = 1 ; Pi{E,a,M) =E-1 


P 2 {E, a,M) = E^ -6E + 5- 8aM 


PsiE, a, M) = E^ - + (59 _ 56aM -f 48o)P -H (120aM - 48a - 45) 

Pa{E, a, M)=E^ - 28E^ + (254 -h 288a - n6aM)E^ - (812 -f 2112a - 1568aM)P -h 
-^(585 -f 1824a - 2160aM - 1920a2M -h 960a‘^M‘^) 

P^iE, a, M)=E^ - AbE'^ + 10(73 -f 96a - 40aM)E^ - (5130 -f 17088a - 7920aM)E^ 
-^(14389 -h 77376a -f 32256a2 - 42032aM - 50304a^M -h 13504a^M2)P - 
(9945 -f 61248a -h 32256a^ - 46800aM - 124032a^M -h 43200a^M^)(3.5) 

The energy eigenvalues are the roots of the polynomial equation such that the 
1 0 D[E) vanishes 

1 0 0 0 0 ..\ 

5-P 1 0 0 0.. 

48(M-l)a 9-E 1 0 0.. 

0 120 (M- 2 )a 13-P 1 0 .. 

0 0 224(M-3)a 17 - P 1 .. 

0 0 0 ./ 


Hill determinant 11: 


D{E) = 


fl-E 
8Ma 
0 
0 
0 

V 0 







If the infinite matrix is truncated at order N, its determinant Dm{E) is solution 
of the recurrence relation 


Dn+i{E) = (4Af + 1 - E)Dn{E) - 8aN{2N - 1)(M + 1 - N)Dn-i{E) (3.6) 
therefore Dm{E) = {—l)^PN{E,a,M). 

Let us insert the formal expansion Eo{a) = 1 + X]^r , where dn are un¬ 
known variables, into the polynomial Pk{E, a, M) and expand in powers of the 
coupling a 

CXD OO 

Pfe(£; = = ^cW[d,M]aP (3.7) 

n—1 n—0 


The coefficients Cp^^[d, M] depend on up to r = p and vanish if dr are the 
perturbative coefficients of the expansion of the ground energy eigenvalue, and 
k > p. This basic property of the polynomials Pk{E,a, M) allows to translate 
the eqs.(3.4) into recursion relations for the coefficients c^\d,M], which allow 
the exact evaluation of d„ (M) in an automated way. The perturbative expansion 
of other, even and odd, energy eigenvalues may also be performed with small 
changes 0 • 

For instance one obtains for the ground state eigenvalue (to save space, I only 
quote the term of order ): 

Eo{a) = 1 - 2Ma - M{2M - 3)0^ - M{AM^ - 15M -f 12)0^ - M{AOM^ - 

5229 

264M2 -h 516M - 297)0^^ - M(28M'^ - 279M^ + M8M^ -M -k 612)a^ - ... 


,M[14739534473769094441397440311/131072 
1784658246229507614806783961/2048 
6366289590321020410675505193/8192 
25467941727184707982391625/128 
5045855866084769519143395/256 
13767162352069961178387/16 
275387290187286326553/16 
150720427787130918 
972416331565443/2 
359298964716 


15258905673899527021780026567/32768 M + 
507856798436368302225354717/512 M® -f 
114826346450120328738858675/256 -f 
17860681260918574154980695/256 -f 
290745864268625747137893/64 -f 
2148046668981444052251/16 -f 
1798144942712999139 M®® -k 
9879984325256823 M®® -f 
16789635150822 

3534526380M®®] -f ... (3.8) 


I checked the coefficients dn{M) up to n = 11 by performing the the regular 
perturbative expansion for logltp^x)], a rather efficient method for not very large 
order. Simple checks of the coefficients dn{M) for higher orders are provided 
by evaluating them for positive integer values of M, where they reproduce the 
easily obtainable expansions of eqs.(2.7). 
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4 Concluding remarks. 

The properties of anharmonic oscillators perturbed by polynomial potentials 
which correspond to quasi-exactly solvable models are peculiar. The pertur¬ 
bative expansions of the eigenvalues have a finite radius of convergence, which 
evades the general situation. In Sect.3, it was indicated that properties of quasi- 
exactly solvable models may be useful to the study of more general non-quasi- 
exactly solvable models. More specifically, one may evaluate in exact, automated 
way, the perturbative expansions of energy eigenvalues. It would be very inter¬ 
esting to know whether the radius of convergence of these expansions collapses 
to zero, as soon as M differs from a positive integer, or there exist other real 
values of M where such radius is finite. If this were the case, quasi-exactly 
solvable models would have the additional merit of suggesting ways of dealing 
with polynomial perturbations. The answer to this question requires standard 
methods of analysis of coefficients of the perturbative expansions which I hope 
to report in a future work. 

After the present letter was completed, I saw the recent paper by M.Znojil 
, which addresses similar issues, with different techniques and an old letter by 
A.V.Turbiner and A.G.Ushveridze ||^ where a subset of the investigation here 
reported in Sect.2 was performed. 
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